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The TP-AGB phase. 
Lifetimes from C and M star counts in Magellanic Cloud clusters 
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Abstract. Using available data for C and M giants with Mt, i < -3.6 in Magellanic Cloud clusters, we derive limits to the 
lifetimes of the corresponding evolutionary phases, as a function of stellar mass. The C-star phase is found to have a duration 
between 2 and 3 Myr for stars in the mass range from ~ 1.5 to 2.8 M Q . There is also an indication that the peak of C-star lifetime 
shifts to lower masses (from slightly above to slightly below 2 M G ) as we move from LMC to SMC metallicities. The M-giant 
lifetimes also peak at ~ 2 M G in the LMC, with a maximum value of about 4 Myr, whereas in the SMC their lifetimes appear 
much shorter but, actually, they are poorly constrained by the data. These numbers constitute useful constraints to theoretical 
models of the TP-AGB phase. We show that several models in the literature underestimate the duration of the C-star phase at 
LMC metallicities. 



1. Introduction 

Since the work by Frogel et al. (1990, hereafter FMB90), the 
thermally pulsing asymptotic giant branch (TP-AGB) phase 
has been recognized to be an important contributor to the inte- 
grated near-infrared luminosity of young and intermediate-age 
stellar populations. Starting from the observed numbers and lu- 
minosities of C and M-type giants in Magellanic Cloud star 
clusters, FMB90 concluded that up to 40 percent of the bolo- 
metric cluster luminosity comes from stars with M\, Q \ < -3.6, 
which are above the RGB tip and hence belong to the TP-AGB. 
An interesting plot by FMB90 showed this fraction as a func- 
tion of SWB cluster type (Searle et al. 1980), which gives ap- 
proximate estimates of cluster ages. Chariot & Bruzual (1991) 
for the first time used this information to calibrate the amount 
of TP-AGB stars added into models of the spectrophotometric 
evolution of galaxies. This was an important piece of informa- 
tion, previously missing, that has contributed to the success of 
such models in interpreting the observed spectra of galaxies. 

Since then several different approaches have been adopted 
to include the TP-AGB phase into evolutionary population syn- 
thesis models. They often use slightly modified versions of 
the original FMB90 plots to constrain the amount of TP-AGB 
stars (see Mouhcine & Lancon 2002; Maraston 2005, and ref- 
erences therein). Maraston (1998, 2005), for instance, follows 
an empirical approach based on FMB90 data to add the ap- 
propriate "fuel consumption" (equivalent to the contribution of 
these stars to the integrated bolometric light) at LMC metal- 
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licities, and then she uses results from uncalibrated TP-AGB 
evolutionary models from Renzini & Voli (1981) to account 
for the basic dependencies with metallicity. Other authors di- 
rectly include extended sets of TP-AGB evolutionary models, 
for a wide enough range of stellar masses and metallicities, 
into their models. This kind of approach goes from adopting 
a very simplified description of the TP-AGB (e.g. Bressan et 
al. 1994) to adding TP-AGB tracks in which the efficiency 
of the third dredge-up is calibrated to reproduce the C-star 
luminosity functions (CSLF) in the Magellanic Clouds (e.g. 
Marigo & Girardi 2001). Different schemes are adopted by 
Bruzual & Chariot (2003) and Mouhcine & Lancon (2002), 
who adopt TP-AGB tracks that distinguish between the C- and 
M-type phases, but which are not directly calibrated on the 
CSLF: Bruzual & Chariot (2003) make a composition of TP- 
AGB tracks from different sources (namely Vassiliadis & Wood 
1993, Groenewegen & de Jong 1993, and Groenewegen et al. 
1995), but do not check the consistency of these composite 
tracks with the observed properties of resolved AGB popula- 
tions. Mouhcine & Langon (2002) instead compute their own 
TP-AGB tracks adopting modified prescriptions but the same 
values of dredge-up parameters as Groenewegen & de Jong 
(1993); they make several different comparisons with data, 
avoiding however the explicit check with the CSLFs that was 
key to Groenewegen & de Jong (1993). 

Owing to the many interrelated aspects of TP-AGB evolu- 
tion and their complex dependence on metallicity and stellar 
mass, our understanding is that the best approach to be used in 
population synthesis should be the inclusion of TP-AGB tracks 
computed in a self-consistent way and directly calibrated us- 
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ing a set of observables from the LMC and SMC. Once cali- 
brated, these self-consistent sets of TP-AGB tracks present the 
potential advantage of better describing the dependence of all 
TP-AGB properties on the stellar metallicity, and hence they 
should work better for metallicities significantly different from 
the LMC and SMC ones. In this context, the main goal of the 
present paper is to translate FMB90 data into useful quantities 

- the TP-AGB lifetimes as a function of stellar mass - for the 
direct calibration of TP-AGB tracks of LMC ad SMC metallic- 
ities. 

Originally, this work has been motivated by the finding by 
Marigo (2002) that all TP-AGB models computed that far have 
used a very improper prescription for their low-temperature 
opacities. She replaced the opacities for solar-scaled chemi- 
cal mixtures - still now the standard choice in the literature 

- by opacities properly evaluated for the chemical mixtures 
of evolving TP-AGB envelopes. One main evolutionary effect 
of variable molecular opacities is the remarkable reduction of 
effective temperatures for C stars, which then causes an ear- 
lier onset of the superwind regime and a reduction of their 
lifetimes, when compared to models computed at fixed solar- 
scaled opacities. This reduction should not affect the previ- 
ous M-type phase, for which no dramatic change in molecular 
opacities is expected. These findings also cause serious doubts 
about the previously-derived behaviours of TP-AGB evolution- 
ary properties - including lifetimes, termination luminosities, 
effective temperatures, etc - with metallicity. It is then evident 
that new grids of TP-AGB tracks are due, and that the C-star 
lifetimes, once derived from empirical data, could provide im- 
portant constraints to them. 

2. Cluster data 

Magellanic Clouds young clusters are clearly the best objects 
to check TP-AGB lifetimes by means of C- and M-star counts: 
They are populous enough to contain appreciable numbers of 
cool giants, have already been searched for them (FMB90 and 
references therein), and at the same time they have reasonably 
well-known distances, ages and metallicities. In order to relate 
the observed numbers of C and M stars to their lifetimes, we 
also need a measure of each cluster's size. The total masses 
are very uncertain even for the best studied clusters, and are 
conditioned by the large mass fraction locked up in low-mass 
dwarfs. In practice, a better measure of cluster's size comes 
from its integrated luminosities. In this paper, we deal with 
the integrated V-band luminosity, Ly, for a series of reasons: 
First the integrated V magnitudes are known for all clusters in 
FMB90's catalog; second Ly smoothly declines with cluster 
age and it is expected to be little sensitive to errors in the clus- 
ter metallicity (cf. Girardi 2000); and third Ly samples stars in 
well-populated evolutionary stages (namely close to the main- 
sequence turn-off and core-He burning, see Chariot & Bruzual 
1991; Girardi & Bica 1993) and hence it is little affected by 
stochastic cluster-to-cluster variations in their number of stars. 
In comparison, the integrated luminosities in red and near-IR 
passbands, like / and K, although available from wide-area sur- 
veys such as 2MASS and DENIS, are too sensitive to the stars 
in the upper part of the RGB and AGB. They present signif- 



icant non-monotonic behaviours with both age and metallic- 
ity, and stochastic cluster-to-cluster variations (Girardi 2002), 
which we prefer to avoid. 

Therefore, a good starting point to derive lifetimes is to use 
the observed total number of C- and M-stars in a cluster, Nc and 
Nm, divided by its integrated V-band luminosity, Ly. These ob- 
served quantities are directly proportional to the M and C-type 
lifetimes. For the M stars, we limit the comparison to the AGB 
stars above the RGB-tip, i.e. those with Mb i < -3.6. Adopting 
LMC and SMC distance moduli of 18.5 and 18.9 mag, this limit 
correspond to entries of mb i < 14.9 and m^ \ < 15.3, for LMC 
and SMC stars respectively, listed in the table 1 of FMB90. 

Table ^ summarises the cluster data available to our pur- 
poses. We have considered all clusters in FMB90, excluding 
the very young ones (i.e. those with t ^ 10 s yr) and a few SMC 
clusters for which we did not find age determinations based on 
the main sequence turn-off photometry. For each cluster, the 
entries in the table correspond to: 

- The age t as derived from main-sequence turn-off photom- 
etry, if available. It is taken from Girardi et al. (1995) for 
LMC clusters, and Da Costa & Hatzidimitriou (1998) and 
Mighell et al. (1998) for SMC clusters. 

- The age-parameter S, as determined and calibrated by 
Girardi et al. (1995) for LMC clusters only. For the most 
populous young and intermediate-age LMC objects (in- 
cluding the entries in Tabled), the relation log(f/yr) = 
6.227 + 0.0733 S gives the turn-off age t with an error of 
about 0.15 dex in log?. For SMC clusters, we derive the 
S parameters directly from the turn-off age using the same 
relation as for LMC clusters. 

- Recent determination of [Fe/H] when available from ei- 
ther spectroscopic data, or from methods directly calibrated 
with spectroscopy; this [Fe/H] list is certainly very hetero- 
geneous and likely incomplete. 

- The C-star counts, and M-star counts above M\, Q \ < -3.6, 
Nc and A^m- They include the spectroscopically confirmed 
C- and M-stars considered to be cluster members - i.e. lo- 
cated within a circle of diameter "close to or somewhat 
larger than 1 arcmin" around each cluster - in FMB90's ta- 
ble 1, and the few objects with dubious membership ("Y?" 
in his column 4) or spectral classification ("C?" or "M?" in 
his column 3). Clusters without such stars are also included 
in the table. 

- The cluster integrated V-band magnitude from Bica et al. 
(1996) for the LMC, and van den Bergh (1981) for the 
SMC. This can be easily converted into the integrated lumi- 
nosity in solar units, after assuming the Sun has a V-band 
magnitude of My Q = 4.847, and apparent distance moduli 
of 18.6 and 19.0 mag, for the LMC and SMC respectively. 

- The surface density of field C-stars, <x c in units of stars 
per deg 2 , at the cluster position, as derived by Blanco & 
McCarthy (1983). We give just the approximate position of 
each cluster in these isopleth maps (their figures 2 and 3). 

For the SMC clusters NGC 419 and NGC 152, FMB90 cat- 
alogue contains many member stars without spectral classifi- 
cation but bright enough (mboi < 15.3) to be in the TP-AGB 
phase. For the sake of homogeneity, we do not include these 
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objects. We suspect they correspond to O-rich TP-AGB stars in 
these clusters, which would otherwise be practically missing. 
This fact warn us that the data for the O-rich TP-AGB stars 
above M\, ] = -3.6 may be severely incomplete for the SMC. 
This does not happen for the LMC, where, owing to smaller 
distance (and maybe to the slightly larger mean metallicities 
and smaller T e ff), stars slightly above m\, \ < 14.9 are, as a 
rule, clearly classified by FMB90 as being of early M subtypes 
(M0-M2). 

The reader will also notice the high degree of standard- 
ization we applied to these data, like a single value of dis- 
tance+reddening for all clusters in each galaxy, and the ap- 
proximative ages. They are necessary because, for the moment, 
there seems to be no satisfactorily homogeneous and updated 
compilation of such data. Apart from these problems, it is im- 
mediately evident from Table ^that, in the end, the AGB data 
is very scarce when we consider individual clusters (see also 
Marigo et al. 1996). Just a handful of LMC clusters host more 
than two confirmed C stars, whereas for the SMC there are 
just two of such clusters. Taken cluster-per-cluster, the statis- 



tics provided by these data would be very poor. Therefore, we 
have added together the cluster data in bins of AS =3, which 
corresponds to age bins of A log t = 0.22 and turn-off mass bins 
of AlogM-ro - -0.10. 

The binned data are presented in Table |2 together with 
the expected mean metallicity [Fe/H] for each bin as given by 
Pagel & Tautvaisiene's (1998) age-metallicity relations (AMR) 
for "bursting models" of the LMC and SMC. These AMRs 
fit well the [Fe/H] data of individual LMC and SMC clusters, 
within the observed scatter of about 0.2 dex (lcr) for a given 
age. Finally, Table|5]presents the estimated turn-off mass, Mjo, 
for the middle of each bin, as derived from its mean log t and 
[Fe/H], and Girardi et al. (2000) evolutionary tracks. The mass 
of the more evolved AGB stars should be slightly higher than 
Mjo but still very close to the tabulated values, with maximum 
differences amounting to just ~ 0.1 M Q . 

How large is the contamination by C and M stars belonging 
to the LMC and SMC fields rather than to the clusters ? We look 
for the answer in Blanco & McCarthy's (1983) study of the sur- 
face distribution of C- and M-type giants across the Magellanic 
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Table 2. Final binned data for C and M stars in Magellanic Cloud clusters 
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Fig. 1. The ratio between the number of C and M stars in clusters, and their integrated V-band luminosity, Nc/Lv and Nm/Lv 
(bottom and top panels, respectively), as a function of turn-off mass. The data are shown separately for the LMC (left panel) and 
SMC (right). 



Clouds. Their figures 2 and 3 show isopleths of field C-star den- 
sities, crc, over both galaxies; the last column of Table^reports 
the location of FMB90 clusters in these isopleth maps. The 
highest values of <r c , of the order of 600 deg~ 2 , are found just in 
the central SMC and LMC bar regions. NGC 2058 is the only 
cluster in our sample located in such a high-density field; for 
it, the cr c = 600 deg~ 2 value would translate into an expected 
number of ~ 0.13 field C stars inside the ~ 2 X 10~ 4 deg~ 2 area 
surveyed by FMB90. This is already a very low expectancy 
value. Most of the clusters in Tabled however, are in outer 
LMC and SMC regions with <r c values well below 300 deg~ 2 . 
Adding all the individual cr c values multiplied by each cluster 
area, the total number of contaminating C stars is estimated to 
be comprised between 0.75 and 1 .4. Compared to our total sam- 



ple of 56 C-stars distributed in 3 1 clusters, this contamination 
is small enough to be neglected. Even more importantly, the 
only individual clusters with a significant probability of being 
contaminated by field C stars, i.e. those with ctq > 300 deg~ 2 , 
are young LMC clusters with S < 32, corresponding to turn- 
off masses higher than 3 M . This is exactly the age interval 
with less observed C-stars for which, as we will see later, just 
upper limits to the C-star lifetimes can be derived. Therefore, 
contamination by field C stars seems not to be a problem 1 . 



1 Cioni & Habing (2003), using DENIS data, find a total number 
of 7750 C-type stars in the LMC, less than the ~ 11000 expected 
from Blanco & McCarthy's (1983) maps. This would be reassuring 
because the field C-star contamination could be even smaller than here 
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Regarding the M-giants, Blanco & McCarthy (1983) data 
cannot be used to estimate the field contamination since it is 
complete only for spectral types later than M5 (i.e. M5+). Cioni 
& Habing (2003, their table 1), from completely independent 
data, find the C/M0+ ratio for stars above the RGB tip to be 
0.30 and 0.27, for the entire LMC and SMC, respectively. From 
these numbers, we can roughly estimate that the contamina- 
tion by field M-giants is just a few times larger (~ 3.3 and 
-3.7 times, for the LMC and SMC) than the one from C- 
giants. Assuming 3.5 as the overdensity factor for both LMC 
and SMC, we can expect from 2.5 to 5 field M-stars contami- 
nating our sample, which contains 68 M stars in total (65 in the 
LMC and 3 in the SMC). We conclude that the field M-giant 
contamination can be safely neglected as well, at least for the 
LMC. For the SMC, the total number of bona-fide M giants in 
clusters (just 3) is so low that field contamination may indeed 
be an issue; however, as we conclude later in Sect.@J this num- 
ber is also low enough to make the SMC M-star data almost 
useless as a constraint to AGB models. 

The FMB90 data refers only to optically visible stars. 
What about the presence of dust-enshrouded, optically ob- 
scured TP-AGB stars? van Loon et al. (2005) present a com- 
prehensive survey and investigation of such stars in Magellanic 
Cloud Clusters, thus representing the mid-infrared counterpart 
of FMB90 survey. Unfortunately, the two cluster samples are 
quite different, van Loon et al.'s (2005) sample of bright IR 
cluster objects (their tables 2 to 7) would add stars to just three 
of our 31 clusters, namely: 1 AGB C-star in NGC 1783, 1 AGB 
C-star and 1 (post-)AGB C-star to NGC 1978, and 2 AGB C- 
stars for NGC 419. It is clear that dust-obscured objects in clus- 
ters constitute a modest fraction of the total numbers of TP- 
AGB stars. We will come back to this point later. 

The final Nc/Ly and Nm/Lv data are presented in Table [2] 
and Figure^together with upper and lower limits given by the 
68% confidence level interval of a Poisson distribution (i.e. lcr 
for the most populated bins). 

3. Lifetimes as a function of mass 

The Nc/Ly and Nm/Lv quantities in Table[2]could already be 
directly used to constrain theoretical models, since they are 
proportional to the typical C- and M-type lifetimes. Large er- 
ror bars are implied by the low number of available objects for 
some of the age bins. Anyway, the LMC data pose well de- 
fined constraints to the C-type lifetime for masses between 1 .5 
and 2.8 M , whereas at higher/lower masses just upper limits 
are derived. For the SMC, error bars are even larger, but three 
age bins have enough statistics to provide useful constraints to 
N c /L v between 1.2 and 2.4 M Q . The N M /L V data for the SMC 

estimated. However, we remark that the Cioni & Habing (2003) clas- 
sification is based on a photometric criterion that likely misclassifies 
a non-negligible fraction (~ 10 %) of C-type stars as M-type ones, 
therefore there may be no real discrepancy in their estimates when 
compared to Blanco & McCarthy's (1983) one. In this paper, we opt 
to use the Blanco & McCarthy (1983) C-star numbers because they 
are based on a spectroscopic classification which is equivalent to the 
one used by FMB90, and which is considered to be complete for the 
C stars above the RGB-tip. 



is in general of low quality, and it is likely severely affected 
by incompleteness; moreover, at SMC metallicities a signifi- 
cant fraction of the O-rich AGB stars above the RGB-tip may 
be of spectral type earlier than M, and then they may be absent 
from FMB90 catalog. As a consequence, one should better not 
using the Nm/Lv ratios derived in the SMC to constrain AGB 
models. 

The N/Lv values can now be converted directly into stellar 
lifetimes as follows. For a given evolutionary stage j, Nj/Ly 
is related to the the lifetime t, through an age-dependent pro- 
portionality constant, that could be computed by using basic 
population synthesis theory. In this paper, we compute the pro- 
portionality constant in a purely numerical way: we take a set 
of AGB models whose AGB lifetimes t A g B (M[,Z) are exactly 
known for all masses and metallicities. Then, we construct 
isochrones and, by simply integrating the stellar number den- 
sity and V-band luminosity, weighted by the initial mass func- 
tion (IMF), along them 2 , we derive the theoretical isochrone 
Nagb/Lv ratios. The empirical lifetimes of the phase j under 
consideration are then given by the ratio between the observed 
Nj/Lv (Table |3 and the corresponding simulated Nagb/Lv, 
times the tagb(M,Z) lifetime that corresponds to the AGB 
stars in that isochrone. By interpolating between models of sev- 
eral masses and metallicities, we properly take into considera- 
tion the variation with metallicity of the main sequence life- 
times - and hence of the evolutionary rate at which stars leave 
the main sequence - and of the integrated V-band luminosity. 
Both effects play a non negligible role in determining the pro- 
portionality constant between Nj/Ly and tj. 

The final results for the lifetimes of C and M-type giants as 
a function of mass are presented in the last columns of Table|2] 
and in Fig. |2] 

In the present work, we have used the same isochrones 
as in Cioni et al. (2006ab) 3 ; they are based on Girardi et al. 
(2000) tracks for the pre-TP-AGB phases, and completed with 
TP-AGB tracks computed on purpose using Marigo's (2002) 
code. As a matter of fact, the Ly values we have used depends 
very little on the particular set of TP-AGB tracks used in the 
isochrones. We have checked that excluding the TP-AGB from 
Cioni et al. (2006ab) isochrones causes the integrated V-band 
magnitudes to increase by less than 0.1 mag at all ages and 
metallicities relevant to this work; 0.04 mag is the typical value 
for this difference. This means that the possible "systematic" 
errors in the lifetimes, caused by possible errors in our TP-AGB 
tracks, are of the order of just ~ 4 %. This is much smaller than 
the errors caused by the poor statistics in the data (see the 68 % 
confidence-level error bars in Fig.Q, so that we consider this 
problem as being of minor importance. 

More relevant may be the errors caused by uncertainties in 
the AMRs: at a given age and for both the LMC and SMC, 
the mean cluster metallicity may be uncertain at a level of 
~ 0.2 dex, as indicated by a series of papers which reached 
somewhat contrasing results for these relations (see e.g. Pagel 

2 This integration is performed assuming Kroupa's (2001) IMF cor- 
rected for binaries, but the results are quite insensitive to the IMF. 

3 These isochrones are available in 
http : //pleiadi . oapd . inaf . it 
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Fig. 2. The lifetimes of the C- and M-type phases (limited to Mt, i < -3.6 in the case of M stars) as inferred from the cluster 
data, as a function of turn-off mass. The data are shown separately for the LMC (left panel) and SMC (right). 



& Tautvaisiene 1998; Dirsh et al. 2000). Keeping the IMF 
fixed, the integrated My changes with metallicity [Fe/H] at a 
rate AMy/A[Fe/H] which is approximately ~ 0.25 mag/dex 
for ages lower than 1.2 Gyr, and ~ 0.7 otherwise 4 . Therefore, 
errors of ~ 0.2 dex in the AMR would translate in errors of 
~ 15 % in the derived lifetimes. Again, this is still smaller than 
the typical errors caused by the poor statistics; on the other 
hand, it is also evident that our results would benefit from a 
better accessment of the AMRs in both the Magellanic Clouds. 

4. Discussion and conclusions 

Despite the large error bars, Tableland Fig.|2]clearly indicate 
that C-star lifetimes have values of about 2 to 3 Myr, for stars 
in the mass interval from ~ 1.5 to ~ 2.8 M (Fig. [5), and for 
metallicities comprised between the -0.3 and -0.7 dex implied 
by LMC and SMC data. There is also an indication that the 
peak of C-star lifetime shifts to lower masses (from slightly 
above to slightly below 2 M G ) as we move from LMC to SMC 
metallicities. The M-giant lifetimes also peak at about 2 M in 
the LMC, with a maximum value of about 4 Myr. In the SMC 
the M-giant lifetimes appear much shorter, but actually they are 
poorly constrained by present data. 

These lifetimes correspond to the optically-visible TP- 
AGB phase. In their study of IR sources, van Loon et al. (2005) 
find that 1.3-3 M stars spend of the order of 10-20 percent of 
their AGB lifetimes as optically-obscured, bright IR objects, 
with mass losses higher than about 10~ 6 M yr-1. Assuming 
the typical lifetimes above the RGB tip to be 1 Myr, they derive 
l-2x 10 5 yr for the duration of this superwind phase. Moreover, 

4 The change in AM V /A[Fe/H] at 1 .2 Gyr is determined by the pres- 
ence of the RGB and red clump at later ages, whose integrated V-band 
light is more sensitive to [Fe/H] than the one from the main sequence. 



since about 3 x 10 5 yr is the time required for loosing the enve- 
lope masses of the sample stars with the observed (superwind) 
mass loss rates, van Loon et al. estimate that 30-70 percent of 
the mass loss of AGB stars occurs as superwind. If we repeat 
the same reasoning using our own estimates for the lifetimes 
in the TP-AGB phase, which are at least twice larger than the 
value used by van Loon et al. (2005), we find that the superwind 
phase may account for all the mass loss during the TP-AGB 
phase. Of course, this indication is anyway very uncertain be- 
cause derived from small numbers of stars. 

We have also verified that several models in the literature 
present C-star lifetimes significantly shorter than the values we 
find for the LMC, as illustrated in Fig. [3] The comparison is 
shown only for a subset of TP-AGB models which has already 
been used, with various procedures, into population synthesis 
of galaxies. 

These models can be roughly classified into two classes: 
the ones previous to Groenewegen & de Jong (1993) are un- 
calibrated, in the sense that they fail in reproducing the C-star 
luminosities observed in the Magellanic Clouds (as discussed 
for the first time by Iben 1981). Renzini & Voli (1981) TP- 
AGB models belong to this class; as it can be seen in Fig. 
they present too low C-star lifetimes as compared to the data, 
and are systematically shifted towards a much higher range of 
masses. Notice that we plot just a subset of Renzini & Voli's 
a — models (their tables la and If); their models with a > 0, 
which take into account the effect of hot-bottom burning in 
more massive stellar envelopes, would present even lower C- 
star lifetimes for Mjo ^ 4 M , It is clear that these models are 
highly discrepant with the data. 

Groenewegen & de Jong (1993) presented the first cali- 
brated TP-AGB models, in which the poorly known parame- 
ters determining the occurrence and efficiency of third dredge- 
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Renzini & Voli (1981), a=0, Case A, 1=0.02 

Renzini St Voli (1981). a = 0. Case A, Z = 0.004 

- -Groenewegen & de Jong (1994), RE mass loss, AMR 
— Groeneweg en & de Jong (1994), VW mass loss, AMR 
Mouhcine Se Lancon (2002), Z=0,008 

Marigo (2(01), logTJ'« = 6.40, \=0.5, Z = 0.008, K fix. 

Marigo (2(02), logT;™ d = 6.40, \=0.5. Z = 0.008, K var. 




Fig. 3. The lifetimes of the C-star phase in the LMC as inferred 
from the cluster data, as a function of turn-off mass (dots with 
error bars), compared to the predictions of some TP-AGB mod- 
els in the literature. The comparison is made with models com- 
puted for the LMC present metallicity (Z 0.008), whenever 
available. 



up events were tuned so as to reproduce the LMC and SMC 
CSLFs. Indeed, all post-1993 models we plot in Fig. [3]have 
positive C-star lifetimes in the right mass range, with the min- 
imum mass for the present of C-stars being close to 1 .5 M , 
and a peak lifetime located somewhere between 2 and 3 M©, in 
rough agreement with the cluster data. This agreement is just 
expected, since fitting the main features of the CSLF requires 
imposing efficient dredge-up for masses as low as 1 .4 Mq, as 
discussed in detail by Groenewegen & de Jong (1993) and 
Marigo et al. (1999). 

However significant discrepancies appear in the C-star 
lifetimes of calibrated models too. For instance, lifetimes in 
Groenewegen & de Jong (1994) tend to be too short, with max- 
imum values of just 0.6 Myr for models assuming Reimers' 
(1975) mass loss formula (RE in the plot), and 1.8 Myr for 
Vassiliadis & Wood's (1993) one (VW in the plot). The dis- 
crepancies are more evident in the mass range between 1.5 
and 2.5 M Q , where model lifetimes fall to less than 0.5 Myr, 
whereas the ones derived here are closer to 2 Myr. Notice 
that the comparison of Fig. [3] is not strictly correct because 
Groenewegen & de Jong (1994) tracks do not include over- 
shooting, whereas the LMC cluster data have been age-dated 
using overshooting models. We estimate that, in order to be put 
Groenewegen & de Jong models in the same scale as the data, 
their masses should be reduced by a factor of about 20 per cent. 
This would not solve the differences in lifetimes. 

The same discrepancy is shared, to a lower extent, by 
Mouhcine & Lancon's (2002) Z = 0.008 models, which present 
the right peak lifetimes of ~ 2.5 Myr, but again have too low 



lifetimes at masses S> 2 M B . Due to the significantly lower 
metallicities of the old LMC clusters, and to the high metal- 
licity dependence expected for the C-type lifetime, this dis- 
crepancy for Mjo 5 2 M e would likely be reduced if the 
comparison was performed with models of smaller metallic- 
ity (say Z = 0.004). Notice that Mouhcine & Lancon (2002) 
adopt the same dredge-up parameters as Groenewegen & de 
Jong (1993), although many of their model ingredients have 
changed. Therefore, the Mouhcine & Lancon (2002) models 
are not strictly calibrated on the CSLF, although their C-star 
luminosities should not be much far from the observed ones. 

Marigo (2001) models do not present these discrepancies in 
their C-star lifetimes at lower masses, and have a peak lifetime 
of 3.4 Myr, which is well compatible with the data. This could 
be considered as a set of TP-AGB models which comply with 
both CSLF and lifetime constraints. However, we know already 
that these models - as well as all of the previously mentioned 
ones - are wrong for a different reason: they do not consider the 
crucial effect of variable molecular opacities as the chemical 
composition changes along the TP-AGB evolution (see Marigo 
2002). The effect of the variable opacities is illustrated by the 
models labelled "Marigo (2002)" in Fig. |2j it can be noticed 
that assuming the same dredge-up parameters as calibrated by 
Marigo (2001), the change from fixed (solar-scaled) to variable 
opacities (K-fix and K-var cases, respectively) causes a reduc- 
tion of ~ 35 percent in the C-star lifetimes. These would be 
still compatible with the cluster data. However, the new mod- 
els have also lower mean luminosities and hence they fail to 
reproduce the CSLF. In fact, updating of the input physics of 
TP-AGB models is not enough, and a re-calibration of dredge- 
up parameters become necessary in this case. 

In conclusion, we have shown that present data for AGB 
stars in LMC clusters represent useful - and so far neglected 
- constraints to the lifetimes of TP-AGB models. Checking for 
these constraints should be especially important if the TP-AGB 
models are to be used into evolutionary population synthesis. 
In fact, using models with the right luminosities (i.e. calibrated 
with the CSLF) but with too low lifetimes would lead to un- 
derestimate the contribution of TP-AGB stars to the integrated 
light of single-burst stellar populations. 

In a following paper, we will present updated TP-AGB 
models computed for the variable-opacity case, in which the 
lifetimes here derived are adopted, together with the observed 
CSLF in the Magellanic Clouds, in the calibration procedure of 
the main model parameters. 
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